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1.  Introduction: 


Formulae  for  sensitivity  analyses  in  linear  programming  problems  are 

generally  confined  to  parametric  variations  of  the  stipulations  vector 

(right-hand  side)  and  the  functional  coefficients.  There  Is  also  no 

great  trouble  in  effecting  studies  for  possible  variations  in  the 

elements  of  vectors  that  are  not  in  the  basis.  This  is  not  to  say  that 

there  has  been  no  effort  or  interest  in  effecting  such  studies  when 

variations  on  one  or  more  basis  vectors  is  to  be  considered.  See  [2] 

1/ 

and  [8].  There  remains  much  to  be  accomplished,  however,  in  devising 

more  straightforward  and  effective  approaches  which  do  not  depend  on 

special  assumptions  with  respect  to  the  basis  matrices  being  considered. 

In  this  paper  we  present  another  possible  approach.  This  is  done 

2/ 

by  generalizing  a  theorem  which  provides  an  exact  inverse  in  the  form 

of  a  product  modification  to  a  given  inverse  when  the  matrix  has  been 

altered  additively  by  a  matrix  of  a  certain  class.  Since  basis  matrices 

consist  of  linearly  independent  column  vectors,  they  always  have  a  left 
3/ 

inverse.  In  what  follows  we  shall  utilize  this  property  and  thereby 
provide  exact  formulae  for  such  sensitivity  analyses  in  the  general  case 
of  linear  programming  matrices,  including  such  cases  as  networks  and 
distribution  (or  transportation)  model  types. 

2.  Mathematical  Theory: 

We  first  introduce  the  following 

See  also  [3]  where  the  problem  of  optimal  alteration  of  structural 
and  other  coefficients  is  considered. 

2 /  See  [1]  and  [7]. 

J3/  Vide,  e.g.  ,  [4]  for  further  discussion  of  left  inverses. 
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Theorem  1:  (or)  Let  the  mxn  matrix  B  have  a  left 

inverse  B  so  that  B  B  -  1,  where  I  la 
the  nxn  Identity  matrix.  (0)  Let  D 
be  an  mxn  matrix  such  that 

1  1  JL 

B  DB  D  -  pB  D  for  some  real  scalar  p. 

Then  a  left  inverse  of  B  +  o  D,  -  p  *  ^  a, 

a  scalar,  is  B*(I  +t  DB*)  -(I  +t  B*D)B*  , 
where 

(1)  T  -  -  ’  •  ^1+P  fr)’1 


Proof:  We  wish  to  show,  by  suitable  choice  of 

t,  that 


(2) 


-  B*(I  + 


DB*) 


(B  +  a  D) 


Because  of  the  symmetry  of  the  product  of 
the  first  two  terms  this  may  be  rewritten 


(3)  I  -  (I  +tB*D)B*  (B  +oD)  -  (I  +tB*D)  (I  +oB*D)  -I  +(Tfcr)B*D  +toB*DB*D 


By  assumption 

(4)  (B*D)2  -  pB*D 

and  so 

1 

(5)  I  -I  +(t  +a  +  Top)  BTD. 

A  JL 

Thus,  for  Bt(I  +tDBt)  to  be  a  left  inverse 
it  suffices  to  set 

(6)  t  +(J  +  t  o-  p  -0 

or 

(7) 


1+pcr 


■  -  o  (1+pa) 


-1 


Q.E.D 
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In  order  to  have  a  usable  form  for  D  which  has  the  assumed  property 

(4)  regardless  of  the  entries  in  B  we  can  take  D  to  be  a  matrix  with 

only  one  non-zero  column.  Evidently  multiplication  of  any  matrix,  A, 

on  the  right  by  D  results  in  a  matrix  with  at  most  one  non-zero  column— 

1/ 

which  appears  in  the  same  position  as  the  non-zero  column  for  D. 

As  may  then  be  verified,  (4)  holds.  It  should  be  noted, however ,  that 
the  value  of  p  may  vary  with  B^. 

For  such  D  we  have 


(8) 
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In  particular,  if  D  has  only  one  non-zero  element  d^,  say,  in  its  r 
row  and  s^  column,  then  the  right-hand  side  of  (8)  reduces  to 


(9) 


d 

r 


0  • . .  a,  .  •  •  0 
lr 


It  might  also  be  further  noted  that  this  product  form  for  the  new 

1  A 

inverse— i.  e.  ,  (1+ t  BD)B— provides  an  alternative  to  the  iteration 

1/ 

procedure  in  the  modified  simplex  method.  This  may  be  evident  on  taking 

o  *1  and  the  non-zero  column  in  D  to  be  the  entering  column  vector  minus 


1/  Similar  formulae  can,  of  course,  be  developed  for  right  inverses 
and  the  usable  forms  for  D  would  then  consist  of  matrices  with 
at  most  one  non -zero  row. 

2/  And  hence  also  to  the  revised  simplex  method.  See  [2]  and  [5]. 
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Che  departing  column  vector  00  that  then  t  ■  -(1+p) 


-1 


3.  Applications  to  Linear  ProRranmlnR: 

We  now  proceed  to  determine  new  simplex  tableau  elements  as  required 

when  moving  from  B  to  B  +  a  D.  Distinguishing  between  old  and  new 

1/ 

quantities  by  Inserting  carets  over  the  latter, 


(10) 


Yj  ■  (B  +  oD)*Pj  ■  (I  +  tB*D)  B*Pj  -  (I  +  tB*D) 


or,  when  D  Is  of  the  form  (8), 


(ID 


•  yj  *  *.j B#d 


where  d1  ■  (d.  ,...,d  )  is  the  non-zero  column,  s,  in  D,  and 
1  in 

T 

Y,  »(y- . . ,y  y  .).  From  (11)  one  may  further  note  that  If  d 

J  lj  8J 

has  only  one  non-zero  element,  say  d^,  then  B^d  is  df  times  the  r^ 

A 

column  vector  of  B  . 

A 

We  next  obtain,  for  the  most  general  permissible  B  D  , 


(12) 


u)T  -  cj  (B  +ctD)*-  cj  B*  (I  +  tDB*)  -  u)T(I  +  t  DB*) 


When  D  has  form  (8),  this  may  be  simplified  to 

(13)  u)T  -  J  +  t  u>Td(B*)  , 

s 

where  (B1^)  13  the  s^  row  of  B^. 

8 

Similarly , 

(14)  z j  -  o)TPj  -  o)T(Pj  +  t  D  Yj)  -  Zj  +  T  0)TD  Y 


1/  The  notation  conventions  used  here  conform  to  the  ones  in  [2]. 
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ln  general,  while  under  the  assumptions  for  (8) 

(15)  Zj  -  +  t  yBj  U)Td. 

We  may  now  obtain  our  desired  formulae.  First,  the  new  value  of 
the  functional  is 

(16)  z  -  z  +  t  uF  DY 

v  O  O  0 

in  general,  and  under  (8),  this  reduces  to 

(17)  z  -  z  +  t  y  a)  d. 

o  o  so 

For  B  +  a  D  to  be  a  dual  feasible  basis,  it  is  necessary  and 
sufficient  that  z^  >  or,  using  the  above  expression, 

(18)  Zj  -  a(l  +  pa)  1  <uTD  Yj  >  Cj  ,  all  j. 

For  the  form  D  specialized  as  (8),  this  becomes 

(19)  Zj  -  <r(l  +  pa)"1  yflj  U)Td  ^  Cj ,  all  j  . 

Primal  feasibility  on  the  other  hand  corresponds  to 

(20)  Yq  -  a(l  +  pa)"1  B*DYq  >  0 
and  thus,  under  (8), 

(21)  Y  -  a(l  +  pa)  1  y  B^d  >  0. 

o  so  — 


Cf.  (11). 


-6- 


We  may  note  that  the  neceaaary  and  sufficient  conditions  for  primal 
and  dual  feasibility  reduce  to  a  system  of  linear  inequalities  to  be 

i 

satisfied  by  the  scalar  a  when  Y  ,  BrD  and  p  are  regarded  as  fixed. 

These  systems  thereby  define  the  range  of  parametric  variation  for  <r 
when  these  conditions  are  to  be  maintained. 

Here  the  developments  have  proceeded  by  reference  to  left  Inverses. 
Analogous  results  may  be  obtained  for  right  inverses  if  (0)  of  Theorem 

1  1  JL 

Is  replaced  by  the  condition  DBT  D1  pDB  for  some  real  scalar  p 

2 

and  BB  -  I.  Convenient  forms  for  D  would  consist  of  matrices  with 
with  only  one  non-zero  row.  The  form  for  the  perturbed  Inverse  would 
again  be  B^(I  +  t  DB^). 

Further  weakening  of  the  hypotheses  of  Theorem  1  seem  to  be  possible 

only  at  the  expense  of  requiring  various  specializing  assumptions  for 

the  matrices  D  and  B.  For  although  Theorem  1  holds  with  left  Inverse 

replaced  by  the  "triplet  property"  A*AA^  ■  A^,  it  need  not  hold  for 

y  jl 

the  von  Neumann-Rao  generalized  inverse  property.  AATA  -  A  where 
B*is  neither  a  left  nor  a  right  inverse.  The  von  Neumann-Rao  property 
Is  the  one  that  is  essential  for  explicit  representation  of  the  general 
solution  and  consistency  of  the  general  linear  system  Ax  -  b. 


1/  See  [4], 
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